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Abstract. - We prove that for a wide class of quantum spin systems with isotropic Heisenberg 
coupling the energy eigenvalues which belong to a total spin quantum number S have upper 
and lower bounds depending at most quadratically on S. The only assumption adopted is that 
the mean coupling strength of any spin w. r. t. its neighbours is constant for all N spins. The 
coefficients of the bounding parabolas are given in terms of special eigenvalues of the N x N 
coupling matrix which are usually easily evaluated. In addition we show that the bounding 
parabolas, if properly shifted, provide very good approximations of the true boundaries of the 
spectrum. We present numerical examples of frustrated rings, a cube, and an icosahedron. 



Introduction. - The study of finite spin systems is not only interesting in its own right 
but has also important applications for understanding so-called "molecular magnets". The 
synthesis of these systems has undergone rapid progress in recent years building on successes in 
coordination and polyoxometalate chemistry . Each of the identical molecular units can 
contain as few as two and up to several dozen paramagnetic ions. The largest paramagnetic 
molecule synthesized to date, the polyoxometalate {Mo72Fe3o}, contains 30 iron ions of spin 
s = 5/2 ||. It appears that in the majority of these molecules the localized single-particle 
magnetic moments ("spins") couple antiferromagnctically and the spectrum is rather well 
described by the Heisenberg model P pT|] . 

Since the dimension of the Hilbert space for N spins of spin quantum number s, given by 
(2s + 1)^, grows rapidly with N and s, the numerical evaluation of all energy eigenvalues may 
be impossible even if the obvious symmetries of the Hamiltonian are exploited . Hence it 
is indispensable to resort to approximation methods, either analytical or numerical. However, 
it appears also desirable to extend the body of rigorous results on spin systems, which could 
serve as a basis or source of inspiration for the development of approximate models. For 
example, the Berezin-Lieb inequalities p2|, [l3| relating spectral properties of the quantum 
systems to those of their classical counterparts provide a foundation for classical or semi- 
classical treatments of spin systems. Among the Berezin-Lieb inequalities are the following 
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rigorous estimates 

(s + l) 2 Egg* < E min < s 2 E^f lc (1) 

for the ground state energy E m - m < of quite general Heisenberg model systems, where 
-^mfn SIC denotes the minimal energy of the corresponding classical spin system described with 
unit vectors. 

In this article we extend the findings of Lieb and Berezin and prove similar inequalities for 
the extremal values E m - ln (S) and -Emax('S') of the energy levels belonging to a given quantum 
number S of the total spin. It turns out that the S'-resolved eigenvalue spectrum is bounded 
by two parabolas. The proof rests mainly on a linear transformation of the spin observables 
according to the eigenbasis of the coupling matrix J. The bounding parabolas depend only 
on iV, s and three special eigenvalues j, jmin? Jmax 

of J. If -E^fn Slc = Nj m i n , which is satisfied 
for all systems where £mf,i SIC has been determined (see [I4|), our estimate s(s + l)i?^|f^ slc 
for S = improves the lower Berezin-Lieb estimate, eq. (II). Hence we conjecture that also 
for arbitrary S our bounds yield the correct s 2 -term of the asymptotic series of E m i n (S) and 
E min (S) for s -> oo. 

We compare the derived bounds with the actual numerically calculated eigenvalue spec- 
trum for various systems. As expected, the bounds are rather close to the actual spectrum if 
s is large, but there is a considerable gap for small s. For special bipartite systems which we 
call "homogeneous Lieb-Mattis systems" we have derived a closer lower bound even for small 
s. 

In addition, we provide numerical evidence that the bounding parabolas, if properly shifted, 
turn out to be very good approximations of the true boundaries of the spectrum, which is 
especially useful for the approximation of low- lying excitations. 

Definitions and simple results. - We consider spin systems with N spin sites, spin quan- 
tum number s and isotropic Heisenberg coupling between all sites /j, and v with coupling 
constants J MI/ . Let (i — 1,2,3) denote the three components of the spin observable s p 

at site /i and, as usual, 

£ (i) = £^> ( 2 ) 

denote the total spin vector and its various components. All linear operators occuring in this 
context operate on a dim — (2s + 1 ^-dimensional Hilbert space TL. The eigenspace of S 2 
corresponding to the eigenvalue S(S + 1) will be denoted by Hs- 
The Heisenberg Hamilton operator can be written as 

H = J^s^ ■ . (3) 
Throughout this Letter the coupling constants J M „ are assumed to satisfy 

Jfiu Jv^jl ^' 3 = ^ , Jfiu^ (4) 

with j being independent of fi. The latter may be viewed as a kind of weak homogeneity 
assumption. The matrix of coupling constants will be denoted by J. Its trace vanishes due 
to J w = 0. Being symmetrical, it has a complete set of (ordered) eigenvalues ji,... ,]n- 
One of them is the row sum j with 1 = -^(L 1, ■ ■ ■ ,1) as the corresponding eigenvector. 
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Let J' denote the matrix JJ restricted to the subspace orthogonal to 1, and j max (j m in) the 
largest (smallest) eigenvalue of JJ'. Since we made no assumptions on the signs of the J^ u , it 
may happen that some of the numbers j, j m in, jmax coincide. In most cases of interest, the 
N x TV-matrix JJ can be easily diagonalized, either analytically or numerically, in contrast to 
the dim x <iim-dimensional Hamilton matrix. We will denote the a-th normalized eigenvector 

Of J by (ci a , . . . , CNa), i. e. 

^ ' Jp,i/Ci/ a = jaC^a, /I, Q, — 1, . . . , N, (5) 
v 

and 

^ c^c^/j = d a/ 3, a,P=l,...,N, (6) 

where we also allow for the possibility to choose complex eigenvectors. Sums over a = 1, . . . ,N 
excluding ay will be denoted by E' \ where aj denotes the index (within the ordered set of all 
eigenvalues) of the eigenvalue j belonging to the eigenvector 1 . 

For later use we will consider a transformation of the spin observables analogous to the 
transformation onto the eigenbasis of J and define 

Definition 1 T a = ^ c^£^, and Q a =f} a -T a , a = 1, . . . , N. 
The inverse transformation then yields 

En = '^2c lia T a , fi=l,...,N. (7) 



In particular, T aj — S/vN. The following lemma follows directly from the definitions: 
Lemma 1 Ns(s + 1) - E m (|m) 2 = E Q Q» = ±S 2 + El Q* ■ 

Bounding parabolas. - Our main result is formulated in the following theorem: 
Theorem 1 1. The following operator inequality holds: 

+ jminNs(s + 1)<H< J —j^S 2 + j max Ns(s + 1) . (8) 

2. Consequently, if \ if) is an arbitrary normalized vector lying in the subspace Tis, the 
following bounds hold for the expectation value of H: 

] —j^S{S + 1) + j min Ns(s + 1) < (<p\H\<p) < J —j^S{S + 1) + j max Ns(s + 1) . 

(9) 



Proof: We rewrite the Hamiltonian in the following form and conclude 

J_ 

N 

*P P 
4- i. V H . 

N 



# = E J^CvpZI ■ fp = E ^PQP = + E ' 3pQp (10) 

(j,i/a>j3 j3 
> Tff+iminE'^ 



P 

j r?2 i • ( at I i 1 \ 1 n2 1 j ~ jmin ^2 



^S 2 + j min lNs(s + 1) - -S 2 I = J -^^S 2 +j min Ns(s + 1) , 
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using (^), the positivity of Qp, and lemma [!} The other inequality follows analogously. ■ 
As a check of our theorem [I] we consider the mean energy E(S) within the subspace Tis- 
Let Ps denote the projector onto Tis and define E(S) = Tr |iTP,g| /Tr|Ps|. Then we 
obtain the following 

Proposition 1 

J - Jmin 5(5 JmmNs ( s + !) < = _L_ (5 ( S + !J _ ^ + !)) (n) 

< J —jp^s(s + 1) + j ma *iV s ( s + 1) . 

Anticipating that formula ( [Tl] ) for £'(5) is valid, which is proven in [jl5|, this proposition 
of course follows from theorem [j], but we want to give an independent proof of the stated 
inequalities for the interested reader. 

Proof: Since 5 < Ns and 1 < N it follows that 



Moreover, 



and, analogously, 



S(S + 1) < Ns(Ns + 1) < N 2 s(s + 1). (12) 
= Tr {J} =j + J2' 3a <j+(N- l) Jmax , (13) 



> j + (N - l) Jmin . (14) 
The first inequality of the proposition is equivalent to 

((N - - j min ) - Nj) S(S + 1) + N 2 s(s + l)(j + (N- l) Jmin ) < (15) 

or 

- {(N - l)j min + j) (5(5 + 1) - N 2 s(s + 1)) < . (16) 



But ( |16| ) follows directly from (12) and (|14|). The second inequality of the proposition is 

proven in a completely analogous way. I 
In the classical limit s — > oo the bounds of theorem [l] assume the values 

E lowci (S) = Lz2^S 2 +j min Ns 2 , (17) 



E uppci (S) = J ] ™* S 2 + 3m ^Ns 2 . (18) 



It is plausible and can be shown analogously to the proof of theorem [l] |l4| that the corre- 
sponding classical Hamilton function assumes only values between these bounds. The classi- 
cal parabolas (|i"7|), ( ]l8| ) intersect in the point S = Ns,E — jNs 2 which corresponds to the 
(quantum or classical) ferromagnetic ground state. For special systems with nearest neighbour 
coupling J, including the A^-polygon (ring), the cube, and the octahedron, it can be shown 
that the bounds Ei oweI (S) and E uppex (S) are assumed for some classical states Q. 

For a special class of bipartite spin systems we can extend theorem [l] and obtain an 
improved lower bounding parabola. 
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Definition 2 Let the set of spin sites be divided into two disjoint subsets A and B such that 
the coupling constants within A or B are < 0, but > between A and B. Moreover, we assume 
that the partial row sum j + of positive entries of J will be constant for all rows. Hence the 
same holds for j~ = j — j + . Further we assume that J is irreducible, i. e. that the spin system 
cannot be decomposed into unconnected parts. 

Then the spin system will be called an HLM-system ("homogeneous Lieb-Mattis system", 
see p6|~^8[j- Further, we will denote by 8 the (positive) difference between the second smallest 
eigenvalue off and j m in- 

It follows that for HLM-systems |.4| = \B\ = N/2. Examples are even rings or cubes with 
suitable coupling. 

Theorem 2 For FILM systems the following operator inequality holds: 

3 —jf^S 2 + jminNs{s + 1) + S(N - 2)8 < H . (19) 

Proof: Due to the HLM- assumption, J has an eigenvector ^^(lj ■■■ , 1 , — 1, . . . ,—1) with 
eigenvalue j~ — j + , if the spin sites are suitably permuted. By the theorem of Gersgorin 
(c.f. |l9]] , 7.2), j~ — j + will be the smallest eigenvalue j m [ n and non-degenerate (theorem of 
Perron- Frobenius, c.f. (ljj, 9.2). It follows that 

Qmin = jl in • T mln = — (§ A - S B ) 2 = — [2{S 2 A + S 2 B ) - S 2 ) , (20) 

where S A = S M <=.4 §a« an( i §,b = ^ ne l ar S es t eigenvalue of Q m i„ within the 

subspace FLs is assumed for the maximal value of the spin quantum numbers S A = Sb = Ns/2, 
hence 

PsQminPs < ^ (^ + 1) - S(S + 1)) = Ns 2 + 2s - ±S(S + 1). (21) 

Thus we may improve the inequality (|i~2| ) by inserting the following term (suppressing the 
Ps) 



^ {Jot Jm'uijQa ^ & I ^ Qa Qmin J 
a \ a J 



(22) 



= S (ns(s + 1) - jjS{S + 1) - Q min ) (23) 

> 5(N-2)s. (24) 

The last inequality follows from (^l|). This completes the proof. ■ 
It may very well be that a derivation of an improved inequality along this line is also pos- 
sible for other than bipartite systems. Note that we have not yet used symmetry assumptions 
up to weak homogeneity (|J). The key point is to obtain better estimates for Q a than those 

already following from (^|). 

Examples. - In the following figures |l|, ||, and || we present examples for the bounding 
parabolas. The energy spectra of fig. [l] have been calculated for the Hamilton operator of spin 
rings with nearest neighbour and next-nearest neighbour interaction 

(N N \ 

£m • Im+i + a J2 Eh- £m+2 , (25) 




Fig. 1 - Energy spectrum of Heisenberg rings with antiferromagnetic exchange between nearest neigh- 
bours and a competing next-nearest-neighbour interaction. Bounding parabolas according to theorem 
[l] are always depicted by thick solid curves, those according to theorem |^ by dashed curves. The thin 
curves depict the approximating parabolas. 



where the sum is understood mod N. In fig. |lj we display spectra for a ring with N = 6 and 
s = 3/2 for various ratios a of the two coupling constants. 

We also consider the Hamiltonian of a system of spins occupying the vertices of a cube 
as well as of an icosahedron with nearest neighbours interaction of constant strength J, see 
fig- 1 

In all examples the bounding parabolas provide rigorous bounds, but unfortunately not 
very narrow ones, especially for small s. Only in the cases of bipartite systems the lower 
bound is impressively close to the boundary of the true spectrum, even if s is small. 

Approximating parabolas. - From the figures it appears that the extremal values E m - ln (S) 
and -E max (>S') of the exact energy spectrum also lie on approximate parabolas. This is not a 
novel observation. It has been noted on several occasions, that for Heisenberg rings with 
antiferromagnetic nearest-neighbour interaction and with an even number N of sites the set 
of minimal energies E m - m (S) forms a rotational band [^0[|2l[]. Further examples, including 
rings for both even and odd N, tetrahedra, cubes, octahedra, icosahedra, triangular prisms, 
and axially truncated icosahedra, suggest that this "rotational band structure hypothesis" 
holds for general Heisenberg spin systems |2^] . The parabolas of our theorem [j] have a similar- 
curvature but a different shift constant compared with the rotational bands. If they are 
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shifted in such a way, that they meet the highest level of the rotational band, which occurs 
for S = Ns, E = Njs 2 , the resulting approximation is very close to the true boundaries of 
the spectrum, compare the thin curves in figures [l], |^, and [| Moreover, it turns out that the 
shifted parabolas exactly meet the minimal and maximal energies of S — Ns — 1 

E mia (Ns - 1) = jNs 2 + 2( Jmia - j)s (26) 
E maK (Ns - 1) = jNs 2 + 2(j max - j)s . 

Identifying the curvature for the lower parabola with (j — jmin)/N provides a very useful 
approximation for the curvature of the true low-lying rotational band, and thus complements 
the derivations of |2^] which are based on the sublattice structure of the spin array. The 
classical counterparts Ei OWCI (S) and E up p eT (S) of the bounding parabolas are also rather close 
to the true boundaries of the spectrum. 

By using these proposed methods one can predict the detailed shape of the spectrum and 
related properties of molecular magnets without diagonalizing a huge Hamilton matrix. 

The overwhelming numerical evidence for the existence of a rotaional band of minimal 
energies [ p2[ 

E min (S,s)* ^ls(S + l) +E (s) (27) 

for general Heisenberg spin systems cannot be mere coincidence. We believe that the derived 
bounding parabolas provide a first step towards proving this general behaviour. 
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Fig. 3 - Energy spectrum of antiferromagnetically coupled Heisenberg rings with N = 5, s = 7/2 
(l.h.s.; j = 2, j min = -(l + \/5)/2, j max = -(l-VE)/2) and N = 7, s = 2 (r.h.s.; j = 2, j min = -1/3- 
V7/3cos(<j))- x /Jj3sm(<j)) = -1.80194, j max = -l/3+2\/7/3 cos(0) = 1.24698; <j) = 1/3 arctan(3\/3)). 
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